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A search for new heavy particles manifested as resonances in two-jet final states is presented. 
The data were produced in 7 TeV proton-proton collisions by the Large Hadron Collider (LHC) and 
correspond to an integrated luminosity of 315 nb" 1 collected by the ATLAS detector. No resonances 
were observed. Upper limits were set on the product of cross section and signal acceptance for 
excited-quark (q*) production as a function of q* mass. These exclude at the 95% CL the q* mass 
interval 0.30 < m q - < 1.26 TeV, extending the reach of previous experiments. 
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Two-jet (dijet) events in high-energy proton-proton 
(pp) collisions are usually described in the Standard 
Model (SM) by applying quantum chromodynamics 
(QCD) to the scattering of beam-constituent quarks and 
gluons. Several extensions beyond the SM predict new 
heavy particles, accessible at LHC energies, that decay 
into two energetic partons. Such new states may in- 
clude an excited composite quark q* , exemplifying quark 
substructure [3-0|; an axigluon predicted by chiral color 
models a flavour-universal color-octet coloron 
or a color-octet techni-p meson predicted by models of ex- 
tended technicolor and topcolor-assisted technicolor [§- 

m- 

Particularly sensitive to such new objects is the 
dijet invariant mass observable, defined as m 33 = 
^J(E 3 ^ + £>) 2 - (pi'i + p 3 ' 2 ) 2 , where E and p are the 
jet energy and momentum, respectively. Several experi- 
ments have examined m 33 distributions in search of new 



12H17|: recently, 1.13 fb 1 of collision data 



resonances 

at the Fermilab Tevatron collider have excluded the ex- 
istence of excited quarks q* with mass 260 < m q * < 
870 GeV This Letter reports the first search by 

the ATLAS experiment [18J at the LHC for such mas- 
sive particles in pp collisions at a centrc-of-mass energy 
of s/s = 7 TeV, based on a data sample corresponding 
to an integrated luminosity of 315 nb _1 . The analysis 
presented here focused on a search for excited quarks be- 
cause of the accessible predicted cross section 0, Q for 
such particles and the benchmark nature of the model 
that allows limits on acceptance times cross section to 
be set for resonant states with intrinsic widths narrower 
than the experimental resolution. 

The analysis technique consisted of a model- 
independent search for a dijet mass resonance on top of 
a smooth and rapidly falling spectrum and relied on the 
measured m 33 distribution to estimate the background 
level to this new possible signal. In the absence of an 
observed new physics signal, upper limits were deter- 
mined on products of cross section (a) and signal ac- 
ceptance (.4) for several q* test masses for a standard set 
of model parameters. 

The ATLAS detector [lij is a multipurpose particle 



physics apparatus with a forward-backward symmetric 
cylindrical geometry and near 4n coverage in solid an- 
gle [3| . The overall layout of the detector is dominated 
by its four superconducting magnet systems, which com- 
prise a thin solenoid surrounding inner tracking detectors 
and three large toroids with an eightfold azimuthal sym- 
metry. 

The calorimeters, which are surrounded by an exten- 
sive muon system, are of particular importance to this 
analysis. In the pseudorapidity region |?7| < 3.2, high- 
granularity liquid-argon (LAr) electromagnetic sampling 
calorimeters are used. An iron-scintillator tile calorime- 
ter provides hadronic coverage in the range \rj\ < 1.7. The 
end-cap and forward regions, spanning 1.5 < |?7| < 4.9, 
are instrumented with LAr calorimetry for both electro- 
magnetic and hadronic measurements. 

The data sample was collected during stable periods 
of 7 TeV pp collisions using a trigger configuration re- 
quiring the lowest-level hardware-based calorimeter jet 
trigger to satisfy a nominal transverse energy threshold 
of 15 GeV [io.]. This trigger had an efficiency greater 
than 99% for events with at least one jet with transverse 
energy higher than 80 GeV. 

Jets were reconstructed using the anti-fcr jet cluster- 
ing algorithm (2l| with a radius parameter R = 0.6. 
The inputs to this algorithm were clusters of calorime- 
ter cells seeded by cells with energy significantly above 
the measured noise. Jet four- vectors were constructed 
by performing a four- vector sum over these cell clusters, 
treating each as an (E, p) four-vector with zero mass. 
These were corrected for the effects of calorimeter non- 
compensation and inhomogeneities by using transverse- 
momentum (pr) and ^-dependent calibration factors 
based on Monte Carlo (MC) corrections and validated 
with extensive test-beam and collision-data studies [20j, 



22j. The m 3 3 observable was computed without unfold- 



ing jets to hadrons or partons. 

In order to suppress cosmic-ray and beam-related back- 
grounds, events were required to contain at least one 
primary collision vertex, defined by at least five recon- 
structed charged-particle tracks, each with a position, 
when extrapolated to the beamline, of \z\ < 10 cm. 
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Events with at least two jets were retained if the highest 
Pt jet (the "leading" jet) satisfied p 3 ^ > 80 GeV and the 
next-to-leading jet satisfied p? > 30 GeV; this ensured 
that the data sample had high and unbiased trigger and 
jet reconstruction efficiencies. Those events containing a 
poorly measured jet with pt > 15 GeV were vetoed to 
prevent cases where a jet was incorrectly identified as one 
of the two leading jets 23J; this affected the event selec- 
tion by less than 0.5%. The two leading jets were required 
to satisfy several quality criteria [23[ and to lie outside de- 
tector regions where the jet energy was not yet measured 
in an optimal way, such as the interval 1.3 < |^ ct | < 1.8. 
Finally, both jets were required to be in the pseudorapid- 
ity region \rf et \ < 2.5, and their pseudorapidity difference 
was required to satisfy 1 77-^ 1 — if 2 1 < 1.3. These cuts, 
which suppress high-mass SM multijet background, were 
determined by performing an optimization of the poten- 
tial signal from q* decays (using a q* mass of 1 TeV) 
compared with the SM background. There were 132,433 
candidates that satisfied these requirements. 

The final event sample was selected by requiring the 
dijet invariant mass to satisfy > 200 GeV in order to 
eliminate any potential kinematic bias in the m n distri- 
butions from the selection requirements on the jet can- 
didates. There were 37,805 events in this sample, which 
formed the m n distribution shown in Fig. [T] 

MC signal events were generated using the excited- 
quark (qg — > q*) production model Q, 0]. The excited 
quark q* was assumed to have spin 1/2 and quark-like 
couplings, relative to those of the SM SU(2), U(l), and 
SU(3) gauge groups, of / = /' = f a = 1, respectively. 
The compositeness scale (A) was set to th e q * mass. 
Signal events were generated using PYTHIA [24j 6.4.21, 
a leading-order parton-shower MC generator, with the 
modified leading-order MRST2007 [2$ parton distribu- 
tion functions (PDFs) and with the renormalization and 
factorization scales set to the mean pt of the two leading 
jets, pythia was also used to decay the excited quarks 
to all possible SM final states, which were dominantly qg 
but also qW, qZ, and 97. The MC samples were pro- 
duced using the ATLAS MC09 parameter tune ^ 
and a GEANT4-based detector simulation [28j. 

Figure [1] shows the predicted signals for q* masses 
of 500, 800, and 1200 GeV, after all selection cuts. 
The signal acceptance (A), which included reconstruc- 
tion and trigger efficiencies near 100%, was found to 
range from - 31% for m q , = 300 GeV to - 48% for 
m q * — 1.7 TeV [29(. The choice of dijet mass bin- 
ning was motivated by the dijet mass resolution of the 
signal. The predicted experimental width ranged from 
cr m „ /m" ~ 11% at m q , = 300 GeV to cr m „ /m" ~ 7% 
at m q * — 1.7 TeV and was dominated by the detector 
energy resolution. 

The background shape was determined by fitting the 



observed spectrum with the function [16| 

f(x) =pi(l -x) P2 x P3+Pilnx , (1) 



whcr 



>/s/s such that /(l) = and /(0) -> +00, 



and ^{1,2.3,4} are f ree parameters. The x Pilnx factor was 
included to describe the high-ro^ part of the spectrum. 
The function in Eqn. Q] has been shown to fit the m n 
observable well in PYTHIA, HERWIG, and next-to-leading- 
order (NLO) perturbative QCD predictions for pp col- 
lisions at = 1-96 TeV [16[. Studies using pythia 
and the ATLAS GEANT4-based detector simulation were 
performed to demonstrate that the smooth and mono- 
tonic form of Eqn. Q] describes QCD-prcdictcd dijet mass 
distributions in pp collisions at y^s = 7 TeV. There is 
good agreement between the MC prediction and the fit- 
ted parametrization in Eqn. [TJ as evidenced by a x 2 P er 
degree of freedom of 27/22 over the dijet mass range 
200 < m« < 1900 GeV. 

The results of fitting the data with Eqn. [T] are shown 
in Fig. [TJ The presence or absence of detectable 
m r} resonances in this distribution was determined by 
performing several statistical tests of the background- 
only hypothesis. A suite of six tests was employed: 
the BumpHunter (30] |. the Jeffreys divergence (HI, the 
Kolmogorov-Smirnovtest, the likelihood, the Pearson x 2 , 
and the TailHunter statistic '32] . The agreement of the 
data with the background-only hypothesis of a smoothly 
varying and monotonic distribution was determined for 
each statistic by calculating the p- value for the data using 
10 3 pseudo-spectra drawn from Poisson variations seeded 
by the results of the fit of Eqn. [T]to the data. The p- value 
of the background-only hypothesis is defined as the frac- 
tion of pseudo-experiments that result in a value of the 
given statistic greater than the value of the same statis- 
tic found by the fit to the data. The results of all six 
tests were consistent with the conclusion that the fitted 
parametrization described the observed data distribution 
well, with p- values in excess of 51%. These observations 
supported the background-only hypothesis. 

In the absence of any observed discrepancy with the 
zero-signal hypothesis, a Bayesian approach was used 
to set 95% credibility-level (CL) upper limits on a ■ A 
for hypothetical new particles decaying into dijets with 
l?^ *] < 2.5. For each of the test masses (indexed by v) 
corresponding to the excited-quark q* predictions, a like- 
lihood function L v was defined as a product of Poisson 
factors computed for each bin (i) of the m n distribution: 



L u (d\b u ,s)=l[- 



Si{u)} 



-[Ki+siiu)} 



(2) 



where di is the observed number of data events in bin 
i, b V i is the background in bin i obtained as described 
below, and Si(v) is the predicted signal added in bin i by 
the signal template; the latter was normalized to the to- 
tal number of predicted signal events s = JV Siiy). For 
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FIG. 1. The data (D) dijet mass distribution (filled points) 
fitted using a binned background (B) distribution described 
by Eqn. [TJ (histogram). The predicted q* signals 0, E| for 
excited-quark masses of 500, 800, and 1200 GeV are over- 
laid, and the bin-by-bin significance of the data-background 
difference is shown. 



each v, the backgrounds in the bins b V i were evaluated 
from a simultaneous five-parameter fit of the signal and 
background distributions to ensure that the background 
determination would not be biased by the presence of 
any signal. The four background parameters were those 
in Eqn. [I] the fifth parameter consisted of the normaliza- 
tion of the predicted v q* signal template. To avoid ac- 
ceptance bias, the lowest q* test mass used was 300 GeV. 
For every q* mass, Eqn. [2] was computed for a range of 
possible signal yields, s, and the resulting likelihood func- 
tion was multiplied by a flat prior in s to give a posterior 
probability density in s. The 95% probability region was 
then determined by integration of the posterior proba- 
bility distribution. This Bayesian technique was found 
to yield credibility intervals that corresponded well with 
frequentist confidence intervals. This was verified by per- 
forming a series of pseudo-experiments to determine, by 
way of a standard frequentist calculation, the coverage, 
or the fraction of times that the 95% Bayesian credibility 
interval contained the true number of signal events. 

The dominant sources of systematic uncertainty, in de- 
creasing order of importance, were the absolute jet en- 
ergy scale (JES), the background fit parameters, the in- 
tegrated luminosity, and the jet energy resolution (JER). 
The JES uncertainty was quantified as a function of 
and ?7 jot , with values in the range 6 ~ 9% 0, IH 
The jet calibration relied on the MC simulation of the 
response of the ATLAS detector; its uncertainty was con- 
strained by varying the ATLAS simulation and from in 



situ information. The systematic uncertainty on the de- 
termination of the background was taken from the uncer- 
tainty on the parameters resulting from the fit of Eqn. Q] 
to the data sample. The uncertainty on a ■ A due to 
integrated luminosity was estimated to be ±11% (35j . 
The JER uncertainty was treated as uniform in px and 
rf ct with a value of ±14% on the fractional pr resolu- 
tion of each jet [36(. The effects of JES, background 
fit, integrated luminosity, and JER were incorporated 
as nuisance parameters into the likelihood function in 
Eqn. [5] and then marginalized by numerically integrating 
the product of this modified likelihood, the prior in s, 
and the priors corresponding to the nuisance parameters 
to arrive at a modified posterior probability distribution. 
In the course of applying this convolution technique, the 
JER was found to make a negligible contribution to the 
overall systematic uncertainty. 

Figure [2] depicts the resulting 95% CL upper limits on 
<7 • A as a function of the q* resonance mass after incorpo- 
ration of systematic uncertainties. Linear interpolations 
between test masses were used to determine where the 
experimental bound intersected with a theoretical pre- 
diction to yield a lower limit on allowed mass. The cor- 
responding observed 95% CL excited-quark mass exclu- 
sion region was found to be 0.30 < m q * < 1.26 TeV us- 
ing MRST2007 PDFs in the ATLAS default MC09 tune. 
Table U shows the results obtained using CTEQ6L1 (j^ 
and CTEQ5L 38] PDF sets. The variations in the ob- 
served limit associated with the error eigenvectors of 
a CTEQ PDF set were found to be smaller than the 
spread displayed in Table |U The excluded regions were 
~ 30 GeV greater when only statistical uncertainties were 
taken into account. The expected limits corresponding to 
the data sample were computed using an analogous ap- 
proach, but replacing the actual data with pseudo-data 
generated by random fluctuations around the smooth 
function described by fitting the data with Eqn. [1] these 
are shown in Fig. [3J with a resulting expected q* mass 
exclusion region of 0.30 < m q * < 1.06 TeV using 
MRST2007 PDFs. As indicated in Table H the two other 
PDF sets yielded similar results, with expected exclusion 
regions extending to near 1 TeV. An indication of the de- 
pendence of the m q * limits on the theoretical prediction 
for the q* signal was obtained by simultaneously vary- 
ing both the renormalization and factorization scales by 
factors of 0.5 and 2, which was tantamount to modifying 
the predicted cross section by approximately ±20%; this 
changed the observed MRST2007 limit of 1.26 TeV to 
1.32 TeV and 1.22 TeV, respectively. 

In conclusion, a model-independent search for new 
heavy particles manifested as mass resonances in dijet 
final states was conducted using a 315 nb -1 sample of 
7 TeV proton-proton collisions produced by the LHC and 
recorded by the ATLAS detector. No evidence of a res- 
onance structure was found and upper limits at the 95% 
CL were set on the products of cross section and signal 
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TABLE I. The 95% CL lower limits on the allowed q* mass obtained using different PDF sets. 



MC Tune 


PDF Set 


Observed Mass Limit [TeV] 
Stat. Syst. Stat, only 


Expected Mass Limit [TeV] 
Stat. © Syst. 


MC09 [27] 


MRST2007 [25] 


1.26 


1.28 


1.06 


MC09' 2. 


CTEQ6L1 [37] 


1.20 


1.23 


0.99 


PerugiaO [39] 


CTEQ5L [38] 


1.22 


1.25 


1.00 



a The MC09' tune is identical to MC09 except for the PYTHIA [24| parameter PARP(82)= 2.1 and use of the CTEQ6L1 PDF set. 
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FIG. 2. The 95% CL upper limit on a ■ A as a function of 
dijet resonance mass (black filled circles). The black dotted 
curve shows the expected 95% CL upper limit and the light 
and dark yellow shaded bands represent the 68% and 95% 
credibility intervals of the expected limit, respectively. The 
dashed curves represent excited-quark a ■ A predictions for 
different MC tunes, each using a different PDF set. 



acceptance for hypothetical new q* particles decaying to 
dijets. These data exclude at the 95% CL excited-quark 
masses from the lower edge of the search region, 0.30 TeV, 
to 1.26 TeV for a standard set of model parameters and 
using the ATLAS default MC09 tune [27|. This result ex- 
tends the reach of previous experiments and constitutes 
the first exclusion of physics beyond the Standard Model 
by the ATLAS experiment. In the future, such searches 
will be extended to exclude or discover additional hypo- 
thetical particles over greater mass ranges. 
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